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ABSTRACT 

A parametric study on cyclic melting and freezing of an 
encapsulated phase change material (PCM), integrated into a 
solar heat receiver, has been performed. The cyclic nature of the 
present melt/freeze problem is relevant to latent heat thermal 
energy storage (LHTES) systems used to power solar Brayton 
engines in microgravity environments. Specifically, a physical 
and numerical model of the solar heat receiver component of 
NASA Lewis Research Center’s Ground Test Demonstration 
(GTD) project was developed. Multi-conjugate effects such as 
the convective fluid flow of a low-Prandtl-number fluid, coupled 
with thermal conduction in the phase change material, 
containment tube and working fluid conduit were accounted for 
in the model. A single-band thermal radiation model was also 
included to quantify reradiative energy exchange inside the 
receiver and losses through the aperture. The eutectic LiF-CaF 2 
was used as the phase change material (PCM) and a mixture of 
He/Xe was used as the working fluid coolant. A modified 
version of the computer code HOTTube was used to generate 
results in the two-phase regime. Results indicate that parametric 
changes in receiver gas inlet temperature and receiver heat input 
effects higher sensitivity to changes in receiver gas exit 
temperatures. 

NOMENCLATURE 
A solid cross-sectional area 

c specific heat of solid regions 

c p specific heat of working fluid 

f geometric view factor 

h,h* specific enthalpy, heat transfer coefficient 
hjf PCM latent heat of fusion 

. L active tube or cavity length 

m mass flow rate 

M total number of axial nodes along tube 

n outer unit normal 


N total number of tubes in receiver 
P wetted perimeter 

q" heat flux 

Q heat transfer rate 

r radial location 

R thermal resistance 

S* geometric shape factor 

t time 

T temperature 

U velocity vector 

W width of PCM region 

X PCM liquid fraction 

z axial location 

k shell loss function 

6 dimensionless working fluid temperature 

p density 

a Stefan-Boltzmann constant 

to turboaltemator TAC speed 

Subscripts 

abs absorbed 

ap aperture 

avg,max average, maximum 

cav receiver cavity 

ch finned-tube fluid channel 

f working fluid 

in, out working fluid inlet, working fluid outlet 
ref reference value 

Superscripts 

n previous time level 

n+ 1 current time level 
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Fig. 1 Thermodynamic cycle for closed Brayton engine integrated with solar heat receiver. 


INTRODUCTION 

The present paper is an extension of a previous paper (Hall, III 
et al., 1997) by the present authors, in which a code validation 
exercise was undertaken. Numerical results generated by the 
computer code HOTTube, which modeled the thermal 
performance of the solar heat receiver component of the Ground 
Test Demonstration (GTD) project at the NASA Lewis Research 
Center, were compared with available experimental results. 

The intermittent nature of solar energy availability for Earth- 
orbit applications presents a particular challenge for space power 
management schemes during traversal into Earth’s shadow 
(eclipse phase). One alternative to photovoltaics with battery 
storage is solar dynamics with latent heat thermal energy storage 
(LHTES) via solar heat receivers. Solar heat receivers are very 
instrumental components in the production of electric power via 
solar dynamic power systems (SDPSs). In a typical operation, 
the SDPS uses: 1) a concentrator to collect and focus the 

incident energy onto the aperture plane of a central receiver, 2) a 
central receiver to collect and distribute, with minimal losses, the 
reflected energy from the concentrator, 3) working fluid tubes 
aligned along the periphery of the receiver to absorb the 
distributed energy as heat, thus, raising the temperature of the 
working fluid (typically a low-Prandtl -number fluid) flowing 
through the tubes, 4) a turbine to expand the high temperature 
working fluid to produce mechanical work via a rotating shaft, 5) 
a compressor to circulate the working fluid through the working 
fluid tubes, and 6) an alternator to convert mechanical shaft 
motion into electric power. A recuperator is added to increase 
the thermal efficiency of the thermodynamic cycle (typically a 
closed Brayton cycle as depicted in Fig. 1). 

Solar heat receivers employing encapsulated phase change 
materials (PCMs) have the advantage over sensible heat receivers 
of requiring less mass while producing higher energy storage 
densities. This, in turn, makes them ideal candidates for energy 
storage in the space environment where temperatures are 


sufficiently high and PCMs with high latent heats of fusion 
become indispensable. 

Investigations on modeling and testing of solar heat receivers 
and/or their subcomponents employing phase change storage 
have appeared in the literature. Some related papers were 
reported by Strumpf and Coombs (1988), Sedgwick (1988, 
1991), Wichner et al. (1988), Wilson and Flanery (1988), 
Kerslake and Ibrahim (1990), Drake (1990), Kerslake (1991), 
Strumpf et al. (1991, 1994), Scarda (1991), and Namkoong et al. 
(1995). Recently, Shaltens and Mason (1996) reported 
experimental results on the operational and thermal performance 
of the NASA Lewis Research Center’s Solar Dynamic Ground 
Test Demonstration (GTD) project. 

In the present paper, a physical and numerical model is 
developed to study the cyclic behavior of the solar heat receiver 
component of the aforementioned GTD system. A parametric 
analysis is undertaken to delineate sensitivities to controllable 
and measurable parameters such as receiver gas inlet temperature, 
receiver heat input, and TAC speed (via an equivalent working 
fluid mass flow rate). Predicting the changes in thermal 
performance as a result of changes in parameters is critical in 
•determining changes in the so-called thermal state-of-charge 
(SOC) of solar heat receivers. Knowledge of the SOC allows for 
better control strategies relating to power management schemes 
during such operations as peak power demand and emergency 
shutdowns with subsequent restarts. 

DESCRIPTION OF THE GTD SYSTEM 

The Ground Test Demonstration (GTD) project is the world’s 
first full scale demonstration of reliable production of electric 
power via solar dynamics (SD) technology. This 
govemment/industry collaboration is carried out in the NASA 
Lewis Research Center’s large thermal/vacuum facility (tank 6). 
This facility is equipped to provide simulated solar flux in high 
vacuum, similar to that which is encountered in Low-earth orbit 
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(LEO). The primary objectives of this project are to demonstrate, 
using flight prototypical components, that system power 
delivered and system efficiency both fall within design target. 
Moreover, most of the hardware used in the GTD system are 
derived from the SD system designed for the Space Station 
Freedom program. 

The 2 kW e (nominal) GTD system consists of an off-axis solar 
concentrator and solar heat receiver with latent heat thermal 
energy storage (LHTES), both of which are integrated with a 
closed Brayton engine (power conversion unit or PCU). A more 
complete description of the GTD system can be found in the 
paper by Shaltens and Mason (1996). 

GTD SOLAR HEAT RECEIVER MODELING 

Problem Description 

The physical geometry for the present study is depicted in Fig. 
2 and the more detailed single-tube model is shown in Fig. 3. It 
is assumed that each tube in the solar receiver is imparted with 
the same incident solar flux; therefore, only a single tube needs to 
be analyzed, with a subsequent summation over all tubes to 
quantify the total receiver thermal performance. 



GTD Specifications 

The specifications for the solar heat receiver of the GTD 
system used for modeling, in which the present numerical results 
reflect, are as follows: 

• effective cavity diameter = 1 .56 ft (47.55 cm) 

• aperture diameter = 7 in ( 1 7.78 cm) 

• active tube length = 2 ft (60.96 cm) 

• canister outer diameter = 1 .78 in (4.52 cm) 

• hydraulic diameter = 0.045 in (0. 1 143 cm) 

• number of tubes in the receiver = 23 

• number of canisters per tube = 24 

• canister material = superalloy Haynes 188 

• PCM = LiF-CaF 2 

• working fluid = helium/xenon (He/Xe) 

The GTD design is essentially a scale-down version of the Space 
Station Freedom SDPS, which was designed to yield a nominal 
alternator power output of 35 kW c . 


Governing Equations 

The generalized integral equation governing the evolution of 
the enthalpy per unit mass for a region (denoted as k) is given by 

f l t Pkh t dV k = -\ s p k h k U k n k dS k - J s q- n k dS k 0) 

which is coupled to the temperature of the kth region by the 
general equation of state 

T k 

K “ Kef = \ c k (J k ')dT k * a c k {r k - T re j ) (2) 

Tref 


Fig. 2 GTD solar heat receiver design (provided by AlliedSignal 
Aerospace). 
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Fig. 3 Schematic of encapsulated PCM tube configuration with 
annular gas flow. 


for constant c k where each of the k regions is identified as 
follows: k=l (outer canister region); k=2 (PCM region); k=3 
(inner canister region); k=4 (working fluid tube region); k=5 
(working fluid region). In the solid regions and liquid PCM 
region (due to the assumption of no convective motion), 

0 k = 6. 


Discrete Representations 

It should be stated at the outset that the containment canister 
outer and inner regions, along with the working fluid tube region, 
are all considered radially lumped. As a result, the radial index j 
corresponding to the representative temperatures of those regions 
are j=l, j=jmax-l, and j=jmax for the outer canister, inner 
canister, and working fluid tube regions, respectively. The 
remaining PCM region is divided into jmax-3 nodes, for each ith 
axial location along the tube. 

Applying Eq. 1 to a control volume at the ith outer canister 
location (k=l) along the tube results in 
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is written for the present configuration as 
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where i is the discrete index corresponding to the axial direction 

(i=l,2, ,M); j is the discrete index corresponding to the radial 

direction (j=l,2, jmax); n is the discrete index for the 

previous time level and n+1 is the discrete index for the current 
time level For the PCM region (k=2), the discrete equation is 
expressed as 
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which is valid in the region 2 < j < jmax-2. Another important 
consideration is that since the zero reference level for the PCM 
enthalpy per unit mass is in the subcooled regime, and is 
somewhat arbitrary, a “flag” based on temperature is used to 
indicate proximity to the melting point. As a result, an additional 
equation governing the fraction of PCM mass in the liquid phase 
is given by 
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( 9 ) 


where P ch is the channel wetted perimeter, m c h is the channel 
mass flow rate, and S* is a geometric shape factor (estimated to 
be in the range 1 < S* < 1 .2), used to account for the degree of 
departure from triangularity of the finned-tube cross-section. 
Finally, the net rate of axially convected enthalpy, used in Eq. 7, 
is expressed as 

<>°> 


Thermal Radiation Model 

An energy balance on the ith node results in the following 
equation governing the net rate of energy absorbed at the ith 
node: 


QabSi “ fif + Qr, “ Qa Pi - Q shell; ^ ^ 

where for i = 1,2, M 


where the liquid fraction for the ith axial canister location and jth 
node is in the range 0 < Xg < 1. This scheme can be considered a 
hybrid between the enthalpy method and the front tracking 
method, since the interface is not tracked explicitly but an 
additional equation is used. Furthermore, the discrete equations 
governing the evolution of enthalpy per unit mass for the canister 
inner region (k=3) and working fluid tube region (k=4) are 
derived as 
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A M+1 
;=1 

Qap, ~ -^°fi-ap{jt ~~ ^tp] 
and the receiver shell loss per node is given by 

U shell; - NM 


( 12 ) 


(13) 


(14) 


In this equation, the function K(T a vg) is recommended by 
Ensworth et al (1996) to be 



For the working fluid region (k=5), a 1-D, quasi-steady model for 
the axial transport of enthalpy per unit mass is used. A modified 
version of the closed-form solution given by (for an isothermal 
wall condition) 


T s -T f 

0(z) = - — L = e 
T s ~T in 


Ph*z 


( 8 ) 


where b = 0.82, = 1860 R (1033 K), and T avg is the 

instantaneous, spatially-averaged canister outer surface 
temperature, given by 


T avg (t) = 


1 


J Ty(z,t)dz 


0 


(16) 
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Modified HOTTube Code 

The numerical results presented in this present paper were 
generated using a modified version of the computer code 
HOTTube, which is a transient, time-explicit, axisymmetric total 
receiver thermal analysis code. HOTTube was initially developed 
by AlliedSignal Aerospace for the Space Station Freedom solar 
heat receiver. For a more complete description of modified 
HOTTube, including grid size and time step requirements, see the 
paper by Hall, III et al. (1997). 

RESULTS AND DISCUSSION 

The results presented in this section pertain to a baseline 
maximum insolation orbit for the orbital altitude corresponding to 
92 minutes total orbit time with 28 minutes of eclipse. This 
maximum insolation of approximately 1.26 suns (1 sun = 1.37 
kW/m 2 ) is incidentally, the maximum output capability of the 
solar simulator lamps of the GTD system. The corresponding rate 
of energy crossing the aperture plane of the receiver during the 
simulated sun period is approximately 12.6 kW. Another point of 
interest is a newly determined SD-Mir design point (determined 
by AlliedSignal Aerospace), which corresponds to a receiver heat 
input of 15.1 kW and TAC operating speed of 53,000 RPM. It 
should be further pointed out that due to the maximum power 
output capability of the GTD solar simulator lamps, a pseudo-SD- 
Mir design point was effected (in the GTD experiment) by 
increasing the sun period to 77 min during the first two cycles to 
get the equivalent amount of energy into the receiver. 

The nominal SD-Mir design point operating speed of the GTD 
engine’s turboaltemator work-producing shaft is 53,000 RPM, 
which yields a working fluid mass flow rate of approximately 
0.3385 lbn/s (154 g/s). A map of TAC speed versus mass flow 
rate is provided by Shaltens and Mason (1996), from which a 
linear curve fit performed by the present authors is given by 

m = 0.0629 + 5.2*1 O' 6 © 

where to is the TAC speed in RPM and m is the working fluid 
mass flow rate in Ibjs. 

Solar Receiver Energy Balance 

As pointed out in Hall, III et al. (1997), an important 
consideration in the numerical prediction of the solar receiver’s 
thermal performance is the instantaneous overall energy balance. 
This is done to account for every unit of energy crossing the 
aperture plane per unit time. Fig. 4 is a plot of the solar 
receiver’s energy budget for the SD-Mir design point. The 
profile for the rate of energy crossing the aperture plane (labeled 

“ Q in ”) resembles a square pulse function, which illustrates the 
switching nature of consecutive sun periods and eclipse periods. 

If the receiver is truly balanced, then the sum (labeled “ Q Jotal ”) 

of the rate of energy extracted by the gas (labeled “ Q gas ”), rate 
of energy lost by reradiation through the aperture and rate of 


energy lost through the receiver shell (labeled “ Qloss ”)> and rate 

of energy stored inside the receiver (labeled “ Q stor ”) should also 

follow this square pulse function. This is indeed the case, as 
shown in Fig. 4, with a maximum error of less than 3 percent. 
This small error is primarily due to the receiver shell heat loss 
approximation. 


Energy Budget for the SD-Mir Design Point 

Q^.15.1 kW TAC Sp**d=53 t OOO RPM 



Fig. 4 Solar receiver energy balance for the SD-Mir design point 
operating conditions. 


SD-Mir Design Point 



Time Since Start of Initial Sunrise (min) 

Fig. 5 Numerical prediction of maximum canister temperature, 
average canister temperature, and receiver gas exit temperature 
over 9 orbit cycles. 
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Receiver Melt History for SD-Mir Design Point 



Fig. 6 Numerical prediction of the total fraction of PCM in the 
liquid phase over 9 orbit cycles. 

Temperature and Melt Fraction Predictions 
Fig. 5 shows the temporal progression of maximum canister 
outer surface temperature, average canister outer surface 
temperature, and receiver gas exit temperature from startup to 9 
orbit cycles. The receiver heat input is fixed at 15.1 kW and the 
TAC speed varies from zero at startup to 53,000 RPM. The 
receiver gas inlet temperature is fixed at 1508 R (838 K). 

The transient operation for the orbital simulations performed to 
reach the SD-Mir design point is detailed as follows: 

1) the startup mass flow rate is fixed at 0.005 \bjs (2.3 g/s) 
until the maximum canister temperature reaches 1900 R 
(1056 K), after which it is increased to 0.3385 lbn/s (154 
g/s) (this occurs during the 2 nd cycle); faster heatup of the 
receiver is effected as a result; 

2) also during startup, the receiver gas exit temperature is fed 
back to the receiver inlet until the receiver inlet reaches 
1508 R (838 K), which is accomplished during the l sl cycle; 
this is an expedient computational scheme that advances the 
numerical solution to the point where higher fidelity 
modeling can be initiated; 

3) during the first two cycles, the total orbit period is 105 min 
with a sun period of 77 min, after which the total orbit 
period is changed to the baseline 92 min with a sun period of 
64 min; 

From the figure, it is shown that balanced orbit conditions are 
reached on the 4 th orbit. A balanced obit is declared when the 
difference between all calculated temperatures is less than 2 R 
(1.1 K) between successive sunrise and sunset conditions. It is 
further observed that over most of the orbit cycles 

Tmax ' > T oul > T avg 


which is what is typically observed during experimental test runs. 
In addition, perusal of the numerical data reveals that the 
maximum canister, average canister, and receiver gas exit 
temperatures at sunset are approximately 2007 R (1 1 15 K), 1904 
R (1058 K), and 1994 R (1 108 K), respectively, for the balanced 
orbits. These temperatures indicate that during the sun period, a 
local superheat event occurred that caused a sharp rise in receiver 
gas exit temperature significantly beyond the PCM melting point 
of 1873 R (1040 K). This is corroborated by the significant rise 
in receiver gas heat input as shown in Fig. 4, where over the sun 
period, the gas heat input increases from a steady value of 
approximately 7.7 kW to approximately 10.2 kW at sunset. 

The corresponding fraction of receiver PCM in the liquid 
phase is illustrated in Fig. 6. Recall that the melting point of LiF- 
CaF 2 is 1873 R (1040 K) and its latent heat of fusion is 340 
Btu/lb m (789 kJ/kg). For the GTD solar heat receiver, the total 
mass of PCM is 53 lb m (24 kg) which corresponds to a maximum 
latent storage capacity of 18,020 Btu (19 MJ or 5.3 kW hr). It is 
observed that the maximum liquid fraction, which occurs on 
sunset for all representative cycles, reaches a maximum of 
approximately 61 percent during the heatup phase and 59 percent 
during the balanced orbit phases. However, at the end of the 
eclipse phase (sunrise), approximately 12 percent of liquid PCM 
remains in the receiver for the balanced orbit cycles. 


Sensitivity to Receiver Gas Inlet Temperature 

Q.^12.6 kW TAC Sp**d=S3,000 RPM 



Fig. 7 Effect of receiver gas inlet temperature on receiver gas exit 
temperature over 9 orbit cycles for the pseudo- SD-Mir design 
point. 
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Sensitivity to Receiver Gas Inlet Temperature 

Q„*»12.6 kW TAC Sp**d=53,000 RPM 



Sensitivity to Receiver Heat Input 

TAC SpMd>53^00 RPU 



Fig. 8 Effect of receiver gas inlet temperature on average canister 
surface temperature over 9 orbit cycles for the pseudo-SD-Mir 
design point. 


Sensitivity to Receiver Heat Input 

TAC Sp*»d= 53,000 RPM 



Fig. 9 Effect of receiver heat input on receiver gas exit 
temperature over 9 orbit cycles. 


Fig. 10 Effect of receiver heat input on average canister surface 
temperature over 9 orbit cycles. 

Parametric Sensitivity 

Shown in Figs. 7 and 8 are the relative sensitivities of receiver 
gas exit temperature and average canister surface temperature to 
receiver gas inlet temperature. Sensitivities are shown for an 
increase in receiver gas inlet temperature from 1508 R (838 K) to 
1608 R (893 K). For this case, the receiver heat input is fixed at 
12.6 kW and the TAC speed is held at 53,000 RPM, which 
corresponds to a gas mass flow rate of 0.3385 \bjs (154 g/s). It 
appears that the receiver gas exit temperature is more responsive 
to the change in receiver gas inlet temperature, as expected. This 
is primarily due to the smaller thermal interaction time scale for 
the gas than for the diffusion time scale through the solid layers 
of canister material and PCM. In fact, over the balanced orbit 
cycles, the 100 R (56 K) increment in receiver inlet temperature 
effects a 1 1 8 R (66 K) increase in receiver gas exit temperature at 
sunset and a 99 R (55 K) increase in average canister surface 
temperature at sunset. These are increases over that which occurs 
at sunset for a receiver inlet temperature of 1508 R (838 K). The 
corresponding changes for sunrise conditions over the balanced 
orbit cycles result in a 202 R (1 12 K) increase in receiver gas exit 
temperature and a 171 R (95 K) increase in average canister 
surface temperature. Also observed in the figures are the 
qualitative changes in the shapes of the temperature profiles, 
which for the higher gas inlet temperature, results in a A-shape 
over the entire orbit period. This is an indication that the system 
moves away from orbital steady-state conditions. 

The sensitivities to changes in receiver heat input are 
illustrated in Figs. 9 and 10. For these cases, the receiver gas 
inlet temperature is fixed at 1508 R (838 K) and the TAC speed is 
held at 53,000 RPM, which corresponds to a gas mass flow rate 
of 0.3385 \bjs (154 g/s). The receiver heat input is varied from 
the pseudo-SD-Mir design point value of 12.6 kW to the actual 
SD-Mir design point value of 15.1 kW. In the balanced orbit 
regime, the additional 2.5 kW results in an increase in receiver 
gas exit temperature of 122 R (68 K) at sunset and an increase of 


NASA TM-107506 7 



194 R (108 K) at sunrise. Similarly, the additional 2.5 kW 
causes an 84 R (47 K) increase in average canister surface 
temperature at sunset and a 112 R (62 K) increase at sunrise. 
In Fig. 9, the qualitative A-shape is retained for the receiver 
gas exit temperature whereas in Fig. 10, there is a departure 
from the A-shape trend for the average canister surface 
temperature. This departure is attributed to a local superheat 
event that occurred in one of the containment canisters. 

CONCLUSIONS 

A physical and numerical model of the solar heat receiver 
component of NASA Lewis Research Center’s Solar Dynamic 
(SD) Ground Test Demonstration (GTD) project has been 
developed. A parametric study was performed to delineate 
changes in receiver gas exit temperature and average canister 
surface temperature due to changes in receiver gas inlet 
temperature and receiver heat input. Results showed that 
receiver gas exit temperatures were more sensitive to changes 
in receiver gas -inlet temperatures. It was also shown that 
changes in receiver heat input effects a higher sensitivity to 
changes in receiver gas exit temperatures. 
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